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ABSTRACT 

Vibrationally excited molecules with sub-mm rotational transitions are po- 
tentially excellent probes of physical conditions near protostars. This study uses 
observations of the v = 1 and v = 2 ro- vibrational modes of HCN (4-3) to probe 
this environment. The presence or absence and relative strengths of these ro- 
vibrational lines probe the gas excitation mechanism and physical conditions in 
warm, dense material associated with protostellar disks. We present pilot obser- 
vations from the Heinrich Hertz Submillimeter Telescope (HHSMT) and followup 
observations from the Submillimeter Array (SMA). All vibrationally excited HCN 
(4-3) v — 0, v — 1, and v — 2 lines were observed. The existence of the three v = 2 
lines at approximately equal intensity imply collisional excitation with a density 
of greater than (10 10 cm -3 ) and a temperature of > 1000K for the emitting gas. 
This warm, high density material should directly trace structures formed in the 
protostellar envelope and disk environment. Further, the line shapes of the v = 2 
emission may suggest a Keplerian disk. This letter demonstrates the utility of 
this technique which is of particular interest due to the recent inauguration of 
the Atacama Large Millimeter Array (ALMA). 



Subject headings: ISM: molecules — stars: formation — stars: individual (AFGL 
2591) — stars: kinematics and dynamics — stars: protostars — submillimeter: 
stars 
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Introduction 



The vast majority of the molecular lines observed at millimeter wavelengths are asso- 
ciated with the lowest vibrational state (v — 0). However, several ro- vibrational molecular 
transitions have been detected i n the ISM and toward s hot molecular cloud cores (CS, H 2 , 
SiO, HC 3 N, CH 3 CN, and HCN Jziurvs fc Turner! Il986l : Eurnerlll987l : Iftolffs et al.ll201l3 lbh\ 
For vibrational excitation to occur these molecules must be in regions of hot, dense gas 
(perhaps associated with a shock) and/or pumped by infrared radiation. Since these exci- 
tation requirements are expected to be fulfilled only nearby to an embedded source, they 
make vibrationally excited molecules with sub-mm/mm wave rotational transitions poten- 
tially excellent probes of the physical conditions very near embedded infrared sources. This 
is especially tru e for those sources w hich suffer so much extinction that they are unobservable 
in the infrared (jWalker et al.lll994l ). 



Initially, vibr ationally excited HCN was used as a probe of physi cal c ondit i ons for astro- 



chemistry studies ( Morris et al.lll971uWilson et al.lll973l ; lMorrislll975l ). In ll994j . lWalker et al. 
used the fact that the critical density and energy above ground for the v = 1 vibrationally ex- 
cited state of CS (7 — 6) and (10 — 9) was high (n~ 10 8 cm -3 and 2000 K respectively) in order 
to probe protostellar disk structure toward IRAS 16293, pioneering this technique. An exci- 
tation analysis of the vibrational lines suggested emission arising from shocks within a hot, 



dense circumstellar disk perhaps due to self-gravitating disk instabilities, e.g., I Adams et al. 
Jl989h . Since this hot, dense gas could be potentially associated with shocks and/or spiral 



density waves in a protostellar accretion disk, vibrational molecular transitions could be used 
clS db probe of disk structure and planet formation. 

We have chosen to observe the v — 0, v — 1 and v = 2 vibrational modes of HCN (4-3) 
rather than CS. The critical densities of HCN are even higher than CS (> 10 10 cm -3 ), but 
the excitation temperature of the v = 1 ro-vibrational transitions is ~ 1050 K, half that of 
CS. This could make HCN an even better probe of dense material in protostellar accretion 
disks, probing material farther from the protostellar object. For gas at the temperatures 
necessary to collisionally excite these ro-vib rational states, we expect higher J- levels are 
preferentially populated ( jWalker et al.l 119941 ) . Thus, the compromise between atmospheric 
transparency and line brightness makes HCN (4-3) the best choice for use as a probe of 
protostellar structure rather than lower lying rotational transitions. 

The selection rules governing the v — 1 and v = 2 transitions help with observ- 
ing efficiency and interpretation. Emission from the (0, 2 2c , 0) and (0, 2 2d , 0) lines is ex- 
pected to be absent du e to Al = ±1 selection rules if purely radiative excitation dominates 
(jZiurys fc Turner! 1 198 w . If collisional excitation of the v = 1 and v — 2 lines is the dominant 
emission mechanism, we would expect to observe both v = 1 lines at approximately equal 
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intensity, and all three v = 2 lines also at approximately equal intensity. Furthermore, the 
extremely high critical density and temperature of these lines allow observation of emission 
that w ould ordinarily be contaminated by foreground and background emission within the 
beam (IWalker et al.lll994j ). With contamination from the envelope removed, interpretation 
of the morphology and kinematics of protostellar sources should be simplified. 

AFGL 2591 is a testbed of high - mass star formation t hat has unusually high HCN 



abundance (~ 10 6 , iBoonman et al.l ( 120011 )). According to iBenz et al.l ( 120071 ) and refer 



ences therein, it can be also be classified as a high mass protostellar object (HMPO) or 
early "Hot Core." HMPOs are distinctive stages of star formation in which the sources 
are bright in the mid-infrared (10-12/xm) and show a central density concentration. In 
this phase, the protostellar object is still deeply embedded in a collapsing dusty envelope. 
These high-mass objects can effectively heat their envelope by FUV radiation. High en- 
ergy radiation from UV photons or X-rays are unobservable directly, but the effects on 
the molecules within the surrounding environment can be observed in the submillimeter. 
van der Tak et al.l (120061 ) provide interferometric observations of the H^O line at 203 GHz 
showing the existence of a rotating flattened structure, which could be evidence for an ac- 
cretion disk ( CeccareUi et al. 201ol). T he central mass of AFGL 2591 is assumed to between 



10 - 16 M (Marengo et all 12 000) 



van der Tak et al.l 120061 ). and the distance to AFGL 



2591 is 1 ±1.0/0.5 kpc ( Ivan der Tak et al.lll999l ). We adopt a distance of 1 kpc in this Letter 
for all analyses. 

This study makes use of submillimeter interferometric observations from the SMA to 
reduce the effects of beam dilution and increase sensitivity over single dish observations. 
The first goal of this study is to determine the primary excitation mechanism for the vibra- 
tional emission: collisional excitation or IR pumping. If the emission arises from collisional 
excitation, it should trace extremely dense, warm gas of the type expected in accretion disk 
structure. Furthermore, we would expect the emission region to be compact, even at in- 
terferometric resolution if produced by dense gas tracing protostellar disk and/or envelope 
structure. The vastly expanded capabilities of the Atacama Large Millimeter Array (ALMA) 
will further increase spatial resolution and sensitivity to permit imaging of a wide variety of 
sources in these ro-vibrational transitions with greatly reduced integration time when com- 
pared to other observatories. The observations of ALGL2591 presented here demonstrate 
the utility of this technique, which will be even more effective and efficient when applied 
with ALMA. 
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2. Observations 

2.1. HHSMT HCN J=4^3 Observations & Data Analysis 

Observations of the v = and v = (0, l lc , 0) lines for AFGL 2591 (a = 20' 1 29 m 24.90 s , 5 = 
40 11 21.0 ) were made at the HHSMT in February 2000 during priority observing time. The 
receiver was tuned to 354.46033 GHz upper sideband (USB). The average T sys was 528 K. The 
on-source integration time was 48 minutes. A second tuning covering the v = (0, l ld , 0) and 
v — 2 lines was made in December 2002. The receiver was tuned to 356.25571 GHz USB. The 
average T sys for this tuning was 1098 K. The on-source integration time was 93 minutes. The 
beam size was ~22" for both tunings. Pointing accuracy was 2 ' and was checked using planets 
and strong continuum sources several times per observing session. Since the observations 
were taken during priority observing time, t^ghz — 0.08 for all observations. For both 
tunings, we used the MPIfR 345 GHz dual polarization SIS receiver, with 1 GHz bandwidth 
Acousto-Optical Spectrometers providing 1.0 MHz (0.87 km s _1 ) spectral resolution. Main 
beam efficiency of this receiver at 350 GHz is vmr ~ 0-45 for unre s olved sources, computed 



via total power measurements of Mars following iMangum et al.l (119931 ) . Calibration was 
performed using the Hot-Sky-Cold method. The wobble switching observing mode was used 
with a 120" chopper throw, and a 50% duty cycle. We subtracted linear fits to the baselines 
for each spectrum, then co-added the spectra, weighted by RMS noise. The final, co-added 
noise levels in the 354.46033 GHz and 356.25571 GHz tunings were 0.010 K and 0.014 K, 
respectively. Reductions were done using the GILDASS software package. 



2.2. SMA Observations & Data Analysis 

Observations of the vibrationally excited spectral lines of AFGL 2591 (a = 20 h 29 m 24.90 s , 
S = 40°ll'21.0") v = 0, v = (0,l lc ,0), (0, l ld , 0), and v =(0, 2°, 0), (0, 2 2c , 0) and (0, 2 2d , 
0) were observed in a single tuning (354.505476 GHz USB) of the SMA in the subcompact 
configuration on September 8, 2008. Our tuning was specifically chosen to place all six of 
the vibrationally excited HCN lines in the USB. The lower side band (LSB) also included 
CO (3— >2) and H 13 CN, but these lines were not the focus of this study. The brightest calibra- 
tors (mwc349a and 2015+371) closest to the source were used as phase calibrators. Uranus 
served as flux calibrator, and the quasar 3c454.3 was used for passband calibration. The 
precipitable water vapor (PWV) provided in the SMA data headers during observations was 
less than 2.5 mm over the 7.2 hours of observations. The field of view of the SMA primary 



1 http://www.iram.fr/IRAMFR/GILDAS 
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beam is 36" at this frequency. Main beam efficiency of the SMA at 354 GHz is ~0.74 for 
unresolved sources^- 

Image reconstruction was do ne using the Mu ltichannel Image Reconstruction, Image 



Analysis, and Display (MIRIAD) (ISault et al. 1995 ) software package with extensions for the 



SMA using the procedures outlined in IZhaol (120091 ) . The continuum image was constructed 
using molecular line-free "chunks" of the correlator. After molecular line subtraction, the 
equivalent bandwidth for the continuum was 1.3 GHz. In order to generate a continuum 
image, we used the MIRIAD INVERT task with uniform weighting, we generated a "dirty 
image" and performed a hybrid Hogbom/Clark CLEAN using a gain of 0.08 with 35000 
iterations to produce a cleaned image. The number of iterations used for cleaning was deter- 
mined systematically by maximizing the map SNR. The "clean" maps were then RESTORed 
by convolution with a gaussian beam, and were then exported to IDL§| for further analysis. 
The FWHM of the synthesized beam is 3.4" x 2.0" (3400 x 2000 AU) with a position angle 
of 37.4°. 

CLEANed three-dimensional image cubes, generated similarly to the process described 
above, were exported to IDL for spectral line extraction. Continuum emission was subtracted 
in the UV data before the clean and invert steps. The spectral line extraction region is 
defined by a two-dimensional gaussian fit to the CLEANed continuum image. The FWHM 
of the extraction region is 4.0" x 2.6" with a 37.4° rotation angle east of north. No smoothing 
was applied to the v=0, v = (0, l lc , 0), and v = (0,l ld , 0) spectral lines, with a spectral 
resolution of 0.273 km s _1 . The v = (0,2°,0), v = (0,2 2c ,0),w = (0,2 2d ,0) spectra were 
boxcar smoothed by 3. The resulting continuum and HCN (4-3) integrated intensity maps 
and associated RMS noise levels are shown in Figure [TJ 



3. Results 

3.1. HHSMT Spectra 

Spectra of AFGL 2591 observed with the HHSMT are shown in Figure [2j We find bright 
v — HCN (4-3) emission and prominent HCN (4-3) v = (0, l lc , 0) (~ 110 mK) emission. 
However, we find a tentative detection of v = Id line (integrated intensity SNR of 7.6, 1.5a) 
to be 5 x weaker than expected (i.e. equal in intensity to the v — lc line). We find no v = 2 
lines in the spectrum. Line parameters and upper limits are given in Table 1. 



smal. sma.hawaii.edu/beamcalc. html 

3 IDL, ENVI, and Image Access Solutions (IAS) are trademarks of Exelis Corporation 
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5.6 " 

3.7" x 2.0" 



Fig. 1. — Integrated intensity contours of HCN emission overlaid on the continuum image in 
greyscale. The peak flux in the continuum is 1.42 Jy beam -1 with an rms noise of 3.4 mJy 
beam -1 . The spectra were extracted from the white dashed ellipse shown in the upper left 
panel. The rms map noise is 5.44, 0.92, 0.89, 0.81, 0.82, and 0.75 K km s -1 from top 
left to bottom right. The top right of each panel is the SNR of the integrated intensity 
maps, defined as (half the maximum value from a gaussian fit to the integrated intensity 
(HM))/(the background noise in the integrated intensity image). The 3.4" x 2.0" synthesized 
beam is the solid ellipse. The red contour levels of the top three panels begin and are spaced 
at 5cr intervals of the rms map noise. The red contour levels of the bottom three panels 
begin and are spaced at 3a intervals of the rms map noise. The lower right hand corner of 
each panel shows the size of the FWHM fit to the integrated intensity, also shown as a blue 
contour. 




3.2. SMA Data 



The v = integrated intensity map in panel 1 of Figure [T] shows evid ence for an 
extended source towards the northwest of the image, consistent with results from lBenz et al. 
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Fig. 2.— Spectra of AFGL 2591 observed with the HHSMT. 



(120071 ) . The v = (0, l lc , 0) and v = (0, l ld , 0) maps show a marginally resolved source. The 
v = (0, 2°, 0), (0, 2 2c , 0) and (0, 2 2d , 0) maps show mostly unresolved sources, with some hint 
of extension in the first 3<r contour level to the southeast in the v = (0, 2°, 0) line. 

Figure [3] shows the HCN (4-3) v = (0,l lc ,0), and v = (0,l ld ,0) spectra, while Figure 
II shows the HCN (4-3) v = (0,2°,0), (0, 2 2c , 0) and (0, 2 2d , 0) spectra. Table □ shows the 
measured line strengths, FWHM, and vi sr of the HCN transitions. The v — transition is the 
brightest with a temperature of ~ 39 K. We observe the v = Id line to be of approximately 
equal brightness to the v — lc line. All three v — 2 lines are of approximately equal 
brightness, with the exception of a single velocity component in the v = 2 2c line. All of the 
HCN emission lines show a v^ sr = —5.3 ± 0.7 km s -1 velocity shift. 

The ratio of the measured brightness temperature of the three v = 2 transitions to the 
two v = 1 transitions can be used to estimate the excitation temperature of the gas in the 
optically thin limit. The ratio of the upper and lower state brightness temperatures are 
related to the excitation energy and excitation temperature for a linear molecule by: 



T, 



MB,! 



MB,/ 



Ai /i/A gi_ e - Ea /T c ^ a 
A/ \vij gi 



X 



where T MB . , are the brightness temperatures (K) of the upper and lower states, A it i are the 
Einstein coefficients, v^i are the upper and lower frequencies, gj / are the upper and lower 
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state degeneracies, Ej/ is the excitation energy, and T ex ,u is the excitation temperature. 
Averaging the excitation temperatures from all six i — 2 to I — 1 line ratios, we obtain an 
estimate of the excitation temperature to be 1050 ± 243 K. 



Benz et all (120071 ) have previously observed AFGL 2591 from 352.6 - 354.6 GHz at 0.6" 
spatial resolution in the extended array configuration. Their observations cover the v — 
and (v — 1) ro- vibrational transitions of HCN. Both groups find a velocity offset of vi sr = 
—5.3 ±0.7 k m s" 1 , consistent with th e systemic velocity of the molecular cloud surrounding 
AFGL 2591 (Ivan der Tak et al.lll999l ). Spatially, we see the sam e extended emission to the 
NW in the v = transition images as Ivan der Tak et al.l (119991 ); Ivan der Tak et al.l (120061 ) 
and Benz et al. J2007h : however, our angular resolution is insufficient to resolve the two 
peaks. They measure a peak brightnesses of 2.3 Jy beam _1 and 3.1 Jy beam -1 , respectively. 
Our measure ments of 24.4 and 6 .6 Jy beam -1 are ~ 11 x and 2.2 x larger than their measured 



values. Since iBenz et al.l (120071 ) use a very different SMA configuration with much higher 
spatial resolution, comparison of the measured line fluxes is not straightforward. The large 
difference in measured flux for the v = HCN line suggests the emission is relatively extended 
and being resolved out in the work of IBenz et al.l (120071 ) . Our larger beam recovers more 
extended flux. The difference between measurements of the v = 1 line fluxes is smaller, 



suggesting marginal resolution of the line flux in IBenz et al.l (120071 ). but to a lesser degree 
than the v = line. 



4. Discussion 

4.1. Origin of Ro- Vibrational HCN Emission 

Observations from the SMA of the v = and both v = 1 lines are ~ 80 x and 90 x 
greater than the HHSMT measurements. For both lines, this is consistent with the ~ 70 x 
reduction in beam dilution provided by the SMA for a source unresolved by both telescopes. 

Given the strong infrared radiation field surrounding AFGL 2591, one might expect 
any detected vibrationally excited HCN emission to be consistent with pumping through 
14 /im absorption. In the SMA data, we observe the all three HCN (4-3) v — 2 lines of 
approximately equal strength. Emission from the v = (0, 2 2c , 0) and v = (0, 2 2d , 0) lines is 
expected to be absent if purely radiative excitation dominates due to the A/ = ±1 selection 
rule. The existence of lines from transitions forbidden to IR pumping at approximately 
the same brightness temperature as a permitted v = (0, 2°, 0) line demonstrates that the 
excitation mechanism for these transitions must be primarily collisional. The small finger of 
emission at the 3a level to the southwest in the v = (0, 2°, 0) line could be evidence of some 



-9 - 




■15 



-10 -5 

Radial Velocity (km/s) 



Fig. 3.— HCN v=0, v=(0,l lc ,0), and v={0,l ld ,0) lines. Velocity resolution is 0.273 km s _1 . 



extended emission resulting from infrared pumping. 

Figure H] shows the v = (0, 2 2c , 0) emission line which is almost twice as strong as the 
other v = 2 lines, with excess emission primarily concentrated in a velocity component at 
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Fig. 4.— HCN v = (0,2°,0), v = {0,2 2c ,0),v = (0,2 M ,0). These data have been boxcar 
smoothed by a factor of n = 3 to a channel width of 0.82 km 



7.5 km s . The source of this excess emission is unclear, and is inconsistent with the 



emission from the other two v — 2 transitions. In a recent SMA study iBenz et al.l (120071 ) 
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note that the HCN (4-3) v — lc line is preferentially emitted very close to the protostar 
(~600 AU) and, unlike other lines observed, its visibility diagram is not well described by 
models including far-UV and X-ray protostellar emission. The excess emi ssion could be 



cause d by maser activity which is seen in other transitions in the source (ITrinidad et al. 



20031 ). although the pumping mechanism is unclear. Another possibility is irradiation by 
FUV/X-ray emission esca ping the central sour ce, enhancing the v = (0, 2 2c , 0) transition in 
the surrounding material (IBruderer et al.l 120091 ) . 



4.2. Double Peaked Emission 



Both the (0, 2°, 0) and (0, 2 2d , 0) spectral lines shown in Figure H] appear to have double 
peaked emission with velocity components at -6.5 km s _1 and -3.5 km s _1 . This could be the 
signature of Keplerian rotation. Construction of position- velocity diagrams are not possible 
because we do not have the spatial resolution necessary to resolve the disk. Instead, we 
measure the separation between the velocity components in the v = (0, 2°, 0) and (0, 2 2d , 0) 
lines using gaussian fitting. This separation is ~ 2.7 km s -1 . Using Kepler's third law, 
assumin g the central mass of AF GL 2591 to be ~ 10 — 16 M Q , a disk inclination angle of 
26 ± 3° (Ivan der Tak et al.ll2006l ). and the velocity centroid separation corresponds to the 
orbital velocity of the outermost disk ed ge, we estimate a disk di ameter of ~ 1000 — 1800 
AU. This result is consistent with that of Ivan der Tak et al.l (120061 ) who estimate the source 
size morphologically using continuum emission measured with a synthesized beam an order 
of magnitude smaller than ours. The possibility also exists that this feature is self -absorption 
i n the emission line, similar to the self absorption in the v = line reported by iBenz et al. 



(120071 ). or result from other kinematic effects like infall. Observation the v = 2 HCN (4-3) 
lines from this source with improved spatial resolution and sensitivity (i.e. with ALMA) will 
be able to clearly prove or disprove the Keplerian disk hypothesis by imaging the kinematic 
structure of the source. 



4.3. Filling Factor 

We have calculated the filling factor of the SMA beam by dividing the brightness tem- 
perature of the spectral lines by the excitation temperature for the spectral feature (see Table 
[I]). This assumes the brightness temperature of the emission region is ~ T ex , true if those 
regions are optically thick. The range of filling factors is as expected, with the v — line the 
largest and the v = 2 lines the smallest. We find the v — emission to approximately fill 
the beam. The filling factors of the v = 1 lines are significantly larger than the v = 2 lines, 
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Table 1: Results (errors) calculated from gaussian fit to the HCN spectral lines. 



Transition 


Frequency 
GHz 


Peak T B SMA 
K 


FWHM 

km s _1 


v lsr 

km s^ 1 


Filling Factor 


l IIIIOIVI ± I 














{U,L ,UJ 


354.460433 


(1.30E-05) 


n i nso (n m r M 

U.1UOU lU.UUUJ 


% 7nnn (f, d.nnnl 


U.OUUU ^O.UUUU ! 




(0,0,0) 


354.505473 


(1.00E-06) 


2.0300 (0.0150) 


5.0000 (0.0700) 


-5.6000 (0.0300) 




(0,2 2c ,0) 


356.135347 


(5.00E-06) 


< 0.0110 








(0,2 2d ,0) 


356.162751 


(5.00E-06) 


< 0.0110 








(0,l ld ,0) 


356.255566 


(1.30E-05) 


0.017 (0.0140) 


8.3000 (4.2000) 


-11.4000 (2.5500) 




(0,2°,0) 


356.301176 


(6.00E-06) 


< 0.0110 








(SMA) 














(0,l lc ,0) 


354.460433 


(1.30E-05) 


10.5240 (0.0305) 


5.1201 (0.0114) 


-5.1294 (0.0073) 


0.0099 (2.8587E-05) 


(0,0,0) 


354.505473 


(1.00E-06) 


38.9784 (0.1511) 


6.6230 (0.0198) 


-4.8467 (0.0126) 


0.9164 (3.5524E-03) 


(0,2 2c ,0) 


356.135347 


(5.00E-06) 


1.9676 (0.0197) 


5.4246 (0.0414) 


-6.4622 (0.0263) 


0.0009 (9.4025E-06) 


(0,2 M ,0) 


356.162751 


(5.00E-06) 


0.8718 (0.0128) 


6.6584 (0.0754) 


-5.6232 (0.0477) 


0.0004 (6.1093E-06) 


(0,l ld ,0) 


356.255566 


(1.30E-05) 


12.2621 (0.0336) 


5.1016 (0.0108) 


-5.2010 (0.0079) 


0.0115 (3.1487E-05) 


(0,2°,0) 


356.301176 


(6.00E-06) 


1.0692 (0.0149) 


7.6095 (0.0847) 


-4.4490 (0.0530) 


0.0005 (7.1861E-06) 



possibly due to shock heating or F UV irradiation of the outer envelope of HCN surrounding 
AFGL 2591 (jBruderer et al.ll2009l ). The filling factors associated with the v — 2 lines are 
approximately equal and very small, suggesting the emission is from compact unresolved 
sources within the beam, highlighting high density, warm structure in the disk and envelope. 



5. Summary 

We have demonstrated the feasibility of a novel technique for the study of protostellar 
accretion disks using vibrationally excited HCN (4-3). We observed AFGL 2591 with the 
HHSMT and the SMA at 350 GHz. Selection rules for the v — 2 HCN transitions allow us to 
determine the excitation mechanism independent of any other observations or assumptions. 
In AFGL 2591, the vibrationally excited HCN emission appears to be collisionally excited. 
The high densities and temperatures required for collisional excitation of these transitions 
imply that the emitting gas should trace structure in the protostellar envelope and/or disk 
surrounding the source. The line shapes of the v = 2 lines suggest the possibility of of an 
unresolved, rotating Keplerian disk. Future followup with ALMA will be able to observe 
AFGL 2591 in the HCN (4-3) v = 2 lines with sufficient sensitivity and spatial resolution to 
clearly detect or rule out the Keplerian disk we suggest. These lines seem to be promising 
probes of protostellar disk and envelope structure. The increased sensitivity and spatial 
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resolution offered by ALMA will allow these tracers to be applied to a variety of protostellar 
sources, and should offer an improved way to isolate and study the very high density, warm 
regions around protostellar objects. 
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